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We draw the three-flavor phase diagram as a function of light- and strange-quark masses for both zero and 
imaginary quark-number chemical potential, using the Polyakov-loop extended Nambu-Jona-Lasinio model 
with an effective four-quark vertex depending on the Polyakov loop. The model prediction is qualitatively 
consistent with 2-1-1 flavor lattice QCD prediction at zero chemical potential and with degenerate three-flavor 
lattice QCD prediction at imaginary chemical potential. 

PACS numbers: 11.30.Rd, 12.40.-y 



Introduction. Determination of the order of QCD phase 
transitions is an important subject not only in hadron physics 
but also in cosmology 1 1]. The chiral and deconfinement tran- 
sitions are widely believed to be crossover at zero chemical 
potential, when physical values are taken for light and strange 
quark masses, m; and rus |0-Hl- However, the order of the 
transitions is sensitive to the number {N f) of flavors and the 
values of m; and mg. A sketch of the three-flavor phase dia- 
gram is plotted in Fig.lTJa) as a function of mi and rris for the 
case of zero chemical potential (p). This sketch, sometimes 
called the Columbia plot, is based on theoretical considera- 
tions and lattice QCD (LQCD) data H-H]. The physical point 
lies near the second-order transition (solid) line. 
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Fig. 1: Sketch of the three-flavor phase diagram in the mi-rris plane. 
Panel (a) shows a sketch for the chiral transition at = 0. The solid 
line denotes the second-order chiral transition line. Panel (b) shows 
a sketch for the RW transition at the end point (T, 9) = (Trw, tt). 
The solid line means the boundary between the first- and second- 
order transition regions. Panel (b) is based on two-flavor lIT^ and 
degenerate three-flavor [17] LQCD results that the RW transition at 
the end point is first order for light and heavy quark masses, but sec- 
ond order for intermediate masses. 

For higher /i, the chiral crossover at the physical point is 
expected to become first order. In this case, there appears a 
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critical end point (CEP) of the first-order transition line, and 
the transition becomes second order on CEP ^t'-'Q] . However, 
clear evidence of the behavior is not shown yet by LQCD be- 
cause of the sign problem at real /i. 

On the contrary, LQCD is feasible [To'-To'l at imaginary 
/i = iOT, where T stands for the temperature and represents 
the dimensionless chemical potential. QCD has a periodicity 
of 27r/3 in 9 called the Roberge-Weiss (RW) periodicity f20l], 
because QCD is invariant under the extended Z3 transforma- 
tion i2lll . At 9 = 7r/3 mod 27r/3, there appears a first-order 
transition at T higher than some temperature Trw [20]. This 
is now called the RW transition ll20ll . On the RW transition 
line starting from the end point (^,T) = j-K.Tfm), a spon- 
taneous C symmetry breaking occurs 1I22I1 . Very recently, 
the order of the C symmetry breaking at the RW end point 
has been analyzed by two-flavor [16^ and degenerate three- 
flavor 1 17] LQCD. For the two cases, the order is first order 
at small and large quark masses, but second order for inter- 
mediate masses. Figure [TJb) is a sketch based on the LQCD 
results for the RW phase transition at the end point . Most of 
the region is unknown at the present stage. 

As an approach complementary to first-principle LQCD, 
we can consider effective models such as the Nambu-Jona- 
Lasinio (NJL) model S [3 and the Polyakov-loop ex- 
tended Nambu-Jona-Lasinio (PNJL) model |j2lll22[ iMjssIl . 
The NJL model can describe the chiral symmetry breaking, 
but not the confinement mechanism. The PNJL model is de- 
signed to make it possible to treat both the mechanisms. The 
effective models have ambiguity in determining their param- 
eters [30|-|32t]. We then take the following strategy. We first 
construct an effective model and determine parameters of the 
model in the regions where LQCD is feasible. Next, we pre- 
dict physical quantities in the regions where LQCD is not fea- 
sible, using the constructed model. 

The original PNJL model cannot reproduce LQCD data 
at imaginary /i quantitatively f3^. This shortcoming of the 
PNJL model seems to be originated in the fact that the cor- 
relation between the chiral condensate a and the Polyakov 
loop is too weak. Therefore, in Ref. fsT], we extended 
the two-flavor PNJL model by introducing the effective four- 
quark vertex depending on ^. This effective vertex includes 
additional mixing effects between a and <!>. The new model 
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is called the entanglement PNJL (EPNJL) model. The two- 
flavor EPNJL model reproduces LQCD data at zero and imag- 
inary fi, particularly on strong correlations between the chiral 
and deconfinement transitions and also on quark-mass depen- 
dence of the order of the RW end point 1 16|. The two-flavor 
EPNJL model reproduces all LQCD data, without chang- 
ing the parameters, at small real /i without and wifli 
strong magnetic field llssll and at finite isospin chemical po- 
tential [34]. 

In this paper, we extend the two-flavor EPNJL model to the 
three-flavor case. Parameters of the three-flavor EPNJL model 
are determined from LQCD data at zero /i and at the RW end 
point. The Columbia pot is drawn for the chiral transition at 
zero /i and for the C symmetry breaking at the RW end point. 

Model setting. We start with the three-flavor PNJL model. 
The Lagrangian density of the model is 
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where aij = (qiqj) and = \/ jP + Mu^ for i, j = u, d, s. 



C =q{i-f^D'' - mo)q + Gs^KqXaqY + [qilbKqf] The dynamical quark mass M,i is defined by 
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where D"" = d" + iA" ^ d" + i5'^gAl\a/2 with the 
gauge coupling g and the Gell-Mann matrices Aq. Three- 
flavor quark fields q — {qu,qd, qs) have current quark masses 
mo = diag(m„, m^, nis). In the interaction part, Gs and Gd 
denote coupling constants of the scalar-type four-quark and 
the Kobayashi-Maskawa-'t Hooft (KMT) determinant inter- 
action [36, 37], respectively, in which the determinant runs in 
the flavor space. The KMT determinant interaction breaks the 
?7a(1) symmetry explicitly. 

In the PNJL model, the gauge field Aj is treated as a homo- 
geneous and static background field ll25ll . The Polyakov loop 
^ and its conjugate <1>* are determined in the Euclidean space 
by 



^* = itrc(Z), 



(2) 



where L — exp(iA4/T) with A4 = iAq. In the Polyakov 
gauge, A4 is diagonal in the color space. The Polyakov po- 
tential U is assumed to be a function of <P and (p* . We take the 
Polyakov potential of Ref. ll27ll : 



a{T) 



b{T) ln(l - 6<Z><P* + 4(<Z>3 + - 3(<Z><I>*)2) 



a(T) = ao + ai[Y 



, (3) 
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Parameters of U are determined to reproduce LQCD data at 
finite T in the pure gauge limit. 

Using the mean field approximation to the quark-quark in- 
teractions in ([T]i, one can get the thermodynamic potential (per 



Mu ^ rrii - 4:Gs(Jii + 2GB(JjjCrkk 
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for 



^ j ^ k. The variables <!>*, <7i{= cr„„ = (7^^) and 
(Ts (= (Tss) are determined by the stationary condition 03211 . 

When /i = i9T, the thermodynamic potential of QCD has 
the RW periodicity ||20|[, i.e. a periodicity of 27r/3 in 9. The 
PNJL thermodynamic potential J? of (|5]) also has this peri- 
odicity, since the potential is invariant under the extended Z3 
transformation il2lll . At 6* = and tt, f2 is C symmetric ll22ll . 
Particularly at 6* = tt, it is spontaneously broken at higher 
T [ 20, ,22] . The order parameter of the spontaneous C sym- 
metry breaking is a 6'-odd quantity such as the imaginary part 
of the modified Polyakov loop ^ = ^e'^ iB . 

An origin of the four-quark vertex Gs is a gluon exchange 
between quarks and its higher-order diagrams. If the gluon 
field A^ has a vacuum expectation value (^o) in its time com- 
ponent, Ai^ is coupled to (Aq) that is related to <P through 
L [38i]. Hence, Gs is changed into an effective vertex Gs(^) 
depending on <P l38t] . Here, the effective vertex Gsi'P) is 
called the entanglement vertex and all interactions including 
Gs(^) are referred to as the entanglement interactions. It 
is expected that dependence of Gs {'P) will be determined 
in the future by the accurate method such as the exact renor- 
malization group method [38-40[. In this paper, however, we 
simply assume the followin g G s {'P) that preserves the chiral 
symmetry, the G symmetry 11221] and the extended Z3 symme- 
try [[2111 : 



Gs{<P) = Gs[l - ai'P'P* - a2{<P^ + <Z>*^)]. 



(7) 



This modification changes the mesonic terms having Gscru 
and the dynamical quark masses Ma in i7. This is the three- 
flavor version of the EPNJL model, and this model has en- 
tanglement interactions in Gs^h and Ma in addition to the 
CO variant derivative included in the original PNJL model. In 
principle, Gd can depend on (p, too. However, we found that 
the (p dependence of Gd yields qualitatively the same effect 
on the phase diagram as that of Gs. As a simple setup, we 
then neglect the (p dependence of Gd • In the present analysis, 
thus, the ^-dependence of Gd is renormalized in that of Gs. 
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In the thermodynamic potential (|5]l, we impose the isospin 
symmetry for the u-d sector (m; = to„ = m^) and take 
the three-dimensional cutoff A for the momentum integra- 
tion 1 32], because this model is nonrenormalizable. Hence, 
the three-flavor PNJL model has five parameters Gs, Gd, mi, 
nis, and A. We use the parameter set of Tablelll lEsI . These pa- 
rameters are fitted to reproduce empirical values of 7r-meson 
mass and decay constant, i^-meson mass and decay constant 
and rf meson mass at vacuum. 

m;(MeV) m^MeV) yl(MeV) GsA^ Gd{0)A^ 
5.5 140.7 602.3 1.835 12.36 

TABLE I: Summary of the parameter set in the NJL sector 1I23I1 . 

Parameters of U are determined to reproduce LQCD data at 
finite T in the pure gauge limit |27]. The original value of Tq 
is 270 MeV, but the deconfinement temperature Tc determined 
by the EPNJL model with this value of Tnis much larger than 
Tc « 160 MeV predicted by full LQCD Therefore, we 

rescale To to 150 MeV so that the EPNJL model can repro- 
duce Tc = 160 MeV. 

The parameters ai and a2 in (Q are so determined as to 
reproduce two results of LQCD at finite T. The first is a result 
of 2h-1 flavor LQCD at ^ = |1] that the chkal ti-ansition 
is crossover at the physical point. The second is a result of 
degenerate three-flavor LQCD at 9 = n that the order 
of the RW end point is first order for small and large quark 
masses but second order for intermediate quark masses. The 
parameter set (ai , a2 ) satisfying these conditions is located 
in the triangle region 

{-1.5ai + 0.3 <a2< -0.86ai + 0.32, > 0}. (8) 

Here, we take ai = 0.25, a2 0.1 as a typical example. 

Results. Figure |2] shows T dependence of light- and 
strange-quark condensates, <ti and <Ts, and the Polyakov loop 

at /X = 0. In the PNJL model of panel (a), ai and as rapidly 
decrease at T w 180 MeV as T increases, after rapidly 
increases at T w 130 MeV as T increases. Thus, the pseu- 
docritical temperature of the chiral crossover is much higher 
than that of the deconfinement crossover The same property 
is also seen in the two-flavor case |30]. In the EPNJL model 
of panel (b), meanwhile, the pseudocritical temperatures of 
the chiral and the deconfinement crossover almost coincide at 
T w 160 MeV. 

Figure [3] shows the order of the chiral transition in the mi- 
nis plane at /i = 0. This figure corresponds to the small rni 
and nis part of FigUfa). The second-order chiral-transition 
line is drawn for three cases, the PNJL result (dotted line) and 
the EPNJL result (solid Une) and LQCD data (+ symbols) 
For each of the three cases, there are the first-order region be- 
low the second-order line and the crossover region above the 
line. The second-order line predicted by the EPNJL model 
is close to that by LQCD data particularly near the physi- 
cal point. Meanwhile, the first-order region predicted by the 
PNJL model is much smaller than that by LQCD data. Thus, 
the EPNJL model yields much better agreement with LQCD 
prediction than the PNJL model. 




0.05 0.1 0.15 0.2 0.25 0.05 0.1 0.15 0.2 0.25 
T (GeV) T (GeV) 



Fig. 2: T dependence of the light- and strange-quark condensates 
and the Polyakov loop at ^ = 0. The quark condensates are normal- 
ized by ai = -0.0142 [GeV^] at T = ^ = 0. Panels (a) and (b) 
represent results of the PNJL and EPNJL models, respectively. 



The deconfinement transitions predicted by the PNJL and 
EPNJL models are crossover in the whole region shown in 
Fig. [3] In the EPNJL model, the crossover deconfinement 
transition almost coincides with the chiral transition, even if 
the chiral transition is crossover 
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Fig. 3: The order of the chiral transition in the mi-iris plane at = 
0. Solid and dotted lines and + symbols represent the second-order 
chiral-transition lines predicted by the PNJL and EPNJL models and 
LQCD 1^), respectively. 

Now we consider the G symmetry breaking at = tt for 
the case of three degenerate flavors {nis — mi)- Figure 2] 
represents the imaginary part of •f' as a function of m; and 
T predicted by the three-flavor EPNJL model. When m; is 
large, the system is close to the pure gauge limit and hence the 
G-symmetry breaking is first order When mi is small, mean- 
while, the system is nearly chirally symmetric and therefore 
the transition is first order In the intermediate mass region, 
the transition is second order The result is consistent with the 
LQCD data 1 17]. 

Figure |5] shows the phase diagram for the G-symmetry 
breaking at the RW end point predicted by the EPNJL model. 
The diagram is plotted as a function of m; and up to A, 
the upper limit for the present model to be applicable. The 
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Fig. 4: The imaginary part of the modified Polyakov loop at S = tt in 
the mi-T plane predicted by the EPNJL model with m — mi — nis. 



two solid lines represent boundaries between the first- and 
second-order transition regions. Below (above) the lower (up- 
per) boundary, the transition is first order The dotted line of 
mi — iris corresponds to the case of Nf = 3. On the dotted 
line, the order is first order for small and large masses but sec- 
ond order for intermediate masses, as expected. At the phys- 
ical point, the order is second order for the present parameter 
set. However, the order can becomes first order at the physical 
point, if we take other parameter sets belonging to the region 
dSj. In the PNJL model, meanwhile, the transition is always 
first order in the entire region of the mi-nis plane. 
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Fig. 5: The order of C symmetry breaking at the RW end point 
predicted by the EPNJL model. The transition is first order below 
(above) the lower (upper) line, while it is second order between the 
two lines. The dotted line stands for a line of mi = ms, that is, the 
case of three degenerate flavors, whereas the x symbol means the 
physical mass. 

In Figs, m and |5] the EPNJL prediction is shown for small 
and large current quark masses niq {q ~ I, s). The appli- 
cability of the NJL-type model to large niq, however, is an 
open question. In fact it was pointed out that nig dependence 
of the chiral-transition temperature is not consistent with the 
corresponding LQCD results Il4li l42ll : as niq increases, the 
chiral-transition temperature goes up sizably in the NJL-type 



model but hardly changes in the LQCD results. In the EPNJL 
model, the chiral-transition temperature almost coincides with 
the deconfinement one that hardly depends on m^, so that the 
EPNJL result is consistent with the LQCD result for the transi- 
tion temperature. It was also pointed out that for large niq the 
pion mass m^r calculated with the NJL-type model is larger 
than the corresponding LQCD result |43]. In the NJL-type 
model the hadron mass calculation is questionable for large 
niq, particularly when the calculated hadron mass is bigger 
than the cutoff A. Therefore, the ENJL predictions shown in 
Fig. |4] and |5] should be regarded as qualitative ones for the 
iTiq > lOOMeV region where the calculated pion mass is big- 
ger than A. However, the fact that there is the second-order 
region at intermediate nig (< lOOMeV) shows that there ex- 
ists a boundary between the first- and second-order regions at 
large niq. In this qualitative sense, the phase diagram of Fig.|5] 
is reasonable for large niq. 

Figure |6] presents the phase diagram in the 6-T plane pre- 
dicted by the PNJL and EPNJL models, where m; and nis 
have physical values. In the PNJL model of panel (a), a first- 
order RW transition (solid) line is connected at the RW end 
point to two first-order deconfinement (dashed) lines. Hence, 
the RW end point is a triple point. In the EPNJL model of 
panel (b), the RW transition is second order at the end point, 
so that there is no first-order deconfinement line connected to 
the first-order RW transition line. For other parameter sets in 
the parameter region (|8), the transition is weak first-order at 
the end point and hence the first-order RW transition line is 
connected at the RW end point to two very-short first-order 
deconfinement lines. 
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Fig. 6: The phase diagram in the 9-T plane predicted by (a) the 
PNJL model and (b) the EPNJL model. Here, physical values of 
mi and m^ are taken. The solid line stands for the first-order RW 
transition line, while the dashed line corresponds to the first-order 
deconfinement line. 



Summary. In summary, we have extended the three-flavor 
PNJL model by introducing an entanglement vertex G's(??). 
The entanglement PNJL (EPNJL) model is consistent with 
2+1 flavor LQCD data for the chiral transition at /i = and 
degenerate three-flavor LQCD data for the RW transition at 
the end point calculated very lately. The three-flavor phase di- 
agram for the RW transition at the end point is first drawn in 
the mi-nis plane by the EPNJL model justified above. 
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